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ABSTRACT 

We present the results of a five year monitoring campaign of the close binary TWA 5 Aab in the TW Hydrae 
association, using speckle and adaptive optics on the W.M. Keck 10 m telescopes. These measurements were 
taken as part of our ongoing monitoring of pre-main sequence (PMS) binaries in an effort to increase the 
number of dynamically determined PMS masses and thereby calibrate the theoretical PMS evolutionary tracks. 
Our observations have allowed us to obtain the first determination of this system's astrometric orbit. We find 
an orbital period of 5.94 ± 0.09 years and a semi-major axis of 0. // 066 ± 0"005. Combining these results 
with a kinematic distance, we calculate a total mass of 0.71 ± 0.14 M (D/44 pc) 3 . for this system. This 
mass measurement, as well as the estimated age of this system, are consistent to within 2a of all theoretical 
models considered. In this analysis, we properly account for correlated uncertainties, and show that while these 
correlations are generally ignored, they increase the formal uncertainties by up to a factor of five and therefore 
are important to incorporate. With only a few more years of observation, this type of measurement will allow 
the theoretical models to be distinguished. 

Subject headings: binaries: visual — stars: pre-main sequence — stars: fundamental parameters — stars: 
individual (TWA 5A) 



1. INTRODUCTION 

Binary stars present a unique laboratory for the study of 
stellar evolution, as their orbital solutions give direct mass es- 
timates. Though the mass of a star is the most fundamen- 
tal parameter in determining its course of evolution, very few 
PMS stars have dynamically determined masses. In these few 
cases, astrometric and spectroscopic studies have lead to PMS 
mass determinations based on the orbital motion of a stellar 
companion or a circumstellar disk relative to the primary star 
(e.g. Ghez et al. 1995, Simon et al. 2000, Steffen et al. 
2001, Woitas et al. 2001, Tamazian et al 2002, Duchene et 
al. 2003). Once measured, these masses can subsequently 
be used to constrain PMS evolutionary models. These mod- 
els have been shown to be systematically discrepant in their 
predictions by up to factors of two in mass and ten in age 
- still, since there are so few well-determined PMS masses, 
little can currently be done to calibrate these tracks. This is 
particularly true in the lowest mass regime, where only three 
systems have total dynamical masses of less than 1M and 
only one single star has a dynamical mass measurement be- 
low 0.5M o (Hillenbrand & White 2004). This study is part 
of an ongoing program to astrometrically determine the or- 
bits of PMS stars and to help constrain these theoretical mass 
tracks. Young stars are particularly important to correctly cal- 
ibrate, as this calibration will aid in the subsequent calibration 
of young brown dwarf and planetary models. 

The TW Hydrae association was originally discovered by 
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Kastner et al. (1997), with only five members confirmed at 
that time. Since its discovery, 25 total members have been 
identified (e.g., Song et al. 2003, Makarov & Fabricius 2001). 
The association has been shown to be quite young (~8-12 
Myr) via lithium abundance tests, space motions, and place- 
ment of its members on the HR diagram (e.g., Zuckerman & 
Song 2004). Additionally, TW Hydrae is quite nearby, with an 
average distance of ~50 pc (e.g., Makarov & Fabricius 2001, 
Mamajek 2005). This makes TW Hydrae one of the closest 
associations of young stars to the Earth and thus is an ideal 
region for studying spatially resolved PMS binaries, as they 
are likely to have orbital periods as short as a few years. 

With these ideas in mind, we began to monitor TWA 5, the 
fifth of the five original members of the TW Hydrae associa- 
tion identified by Kastner et al. (1997). TWA 5 is composed 
of at least three components. The pair that we analyze here is 
TWA 5Aa-Ab, which had a separation of ~0. // 06 when it was 
discovered by Macintosh et al. (2001). TWA 5 A also has a 
brown dwarf companion, TWA 5B, located ~2" away (Webb 
et al. 1999, Lowrance et al. 1999, Neuhauser et al. 2000, 
Mohanty et al. 2003). Finally, TWA 5Aab is suspected to 
contain at least one spectroscopic pair based upon large radial 
velocity variations (Torres et al. 2003, Torres 2005 private 
communication) . 

In this paper, we present the results of six years of speckle 
and adaptive optics (AO) observations of TWA 5 Aab with the 
W.M. Keck 10 m telescopes. In §2, we describe our data re- 
duction techniques, and in §3 we present our orbital solution 
for the system. In §4, we discuss comparisons of our derived 
orbital parameters with mass and age predictions from the- 
oretical PMS tracks and make recommendations for future 
studies of this system. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Speckle Data 

TWA 5 Aab was observed using the Keck I 10 m telescope 
with the facility Near Infrared Camera (NIRC, Matthews & 
Soifer 1994, Matthews et al. 1996) roughly once a year from 
2001 to 2005 (see Table 1 for exact dates). In its high angu- 
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Fig. 1. — The two-dimensional visibilities for all the speckle data are shown in combination with plots of one-dimensional cuts through the visibilities. The 
points represent actual data, while the solid lines show the best fit of Equation 1 to this data. The data from 2002 clearly illustrates that the separation of the 
binary during this epoch was below the diffraction limit. The degeneracy of the separation and flux ratio in this fit (where the first minimum was not reached) 
necessitated the fixing of the flux ratio in this epoch to obtain the correct separation. The results of this analysis can be seen in Figure 3. The slight discrepancy in 
the 2004 fit comes from discrepancy between the calibrator and the object, stemming from large scale changes in seeing on short timescales. However, the large 
error bars on this data point account for this variation and thus they do not weigh heavily in our orbital fit. 



TABLE 1 
TWA 5A Binary Star Parameters 



Date of 


Filter 


Xo 


Separation 


Position Angle 


Flux Ratio 


Speckle 


Source a 


Observation 




(/xm) 


(arcseconds) 


(degrees) 


(Aa/Ab) 


or AO? 




2000 Feb 20 


H 


1.648 


0.0548 ± 0.0005 


205.9 ± 1.0 


1.09 ±0.02 


AO 


2 


2000 Feb 22 


K' 


2.127 


0.054 ± 0.003 


204.2 ± 3.0 


1.11 ±0.07 


AO 


3 




H 


1.648 






1.09 ±0.08 


AO 


3 




J 


1.26 






0.94 ± 0.05 


AO 


3 


2001 May 06 


K 


2.2 


0.0351 ±0.0002 


12.67 ± 1.06 


1.26 ±0.09 


Sp 




2002 May 23 


K 


2.2 


0.013 ± 0.003 


313.66 ±2.99 


1.24±0.08 b 


Sp 




2003 Dec 05 


K 


2.2 


0.0306 ± 0.0004 


227.41 ± 5.49 


1.22 ±0.04 


Sp 




2004 Dec 18 


K 


2.2 


0.0515 ±0.0009 


32.10 ±2.22 


1.39 ±0.09 


Sp 




2005 Feb 16 


Fell 


1.65 


0.053 ± 0.001 


32.59 ± 5.22 


1.29 ±0.18 


AO 




2005 May 27 


K 


2.2 


0.0574 ± 0.0003 


29.68 ± 0.35 


1.23 ±0.04 


Sp 




2005 Dec 12 


Kp 


2.124 


0.0571 ± 0.001 


29.99 ± 2.26 


1.10 ±0.05 


AO 






H 


1.633 


0.0571 ± 0.002 


28.89 ± 0.99 


1.09 ±0.03 


AO 






J 


1.248 


0.0568 ± 0.005 


28.41 ±3.15 


1.05 ±0.13 


AO 




a 1 = This work; 2 


= Macintosh et al. 2001; 2 


= Brandeker et £ 


il. 2003 b Flux ratio fixed in this epoch 









lar resolution mode, NIRC has a pixel scale of 20.44 ± 0.03 
mas/pixel (see Appendix A for details of NIRC's pixel scale 
and orientation). Three to four stacks of 190 images, each 
0.137 seconds long, were obtained through the K band-pass 
filter (X = 2.2 /im, AA = 0.4 /im). They were taken in sta- 
tionary mode, meaning the Keck pupil is fixed with respect to 
the detector during the observations, but the sky rotates. The 
rotation rate is sufficiently slow that it is negligible for indi- 
vidual exposures, but not over the entire stack. Stacks of four 
dark frames were also taken with each object stack. Addition- 
ally, identical stacks of a point source calibrator star (TWA 1 
for 2001-2002 observations, TWA 7 for 2003-2005 observa- 
tions) and an empty portion of sky were obtained immediately 
before or after the target stacks. 

The stacks are processed using image reduction and speckle 
data anaylsis techniques. Specifically, each image is first dark 



and sky subtracted, flat fielded, and bad pixels repaired. The 
images are then individually corrected for a minor optical dis- 
tortion in the NIRC camera (J. Lu et al. in prep). The object, 
calibrator, and sky stacks are then Fourier transformed and 
squared to obtain stacks of power spectra. Next, the calibrator 
stacks and sky stacks are each averaged together (without ro- 
tation). The extraction of the object's intrinsic power spectra 
utilizes the convolution theorem, which gives the relation: 



= 



I-<S> 
<C>-<S>' 



(1) 



where O is the object power spectrum, I is the initial squared 
Fourier transform of the object, <S> is the averaged, squared 
power spectrum of the sky, and <C> is the averaged, squared 
power spectrum of the calibrator. Each power spectrum of the 
object is then rotated so that north is up; these rotated power 
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spectra are combined to obtain an average power spectrum for 
the stack. Finally, those individual stacks were also averaged 
together to produce a final power spectrum. This sequence is 
necessary due to the lack of azimuthal symmetry of the Keck 
I pupil. 

Since a binary star in the image domain is essentially two 
delta functions, the power spectra can be approximated as a 
sinusoid with the functional form: 



p( A = R 2 + l+2Rcos(27rf.f) 
{J) R 2 + \+2R 



(2) 



where R is the flux ratio of the binary star and *f is the vector 
separation of the two stars. This function is fit to our two- 
dimensional power spectra over all spatial frequencies from 
2.68 arcsec -1 to 17.6 arcsec -1 (this procedure is described in 
detail in Ghez et al. 1995). The lower cutoff is imposed to 
avoid spatial frequencies that are corrupted by small changes 
in the atmospheric conditions between the object and the cal- 
ibrator, whereas the upper cutoff is imposed to reject excess 
noise at the highest spatial frequencies. This fitting procedure 
gives a very accurate estimate of the binary star's separation 
and flux ratio for components separated by more than A/2D, 
leaving only a 180° position angle ambiguity. This ambigu- 
ity is resolved by reanalyzing the original images using the 
method of shift and add. The brightest speckle in the speckle 
cloud in each image of the stack is shifted to a common po- 
sition, and then all the images are added together to produce 
a diffration-limited core surrounded by a large, diffuse halo. 
These images allow us to determine the correct orientation 
of the position angles of TWA 5 A. In one epoch (2002), the 
binary separation was less than A/2D, and hence the first min- 
imum of the power spectrum was not measured; without this 
measurement, the separation and the flux ratio are degenerate. 
For this epoch, we fix the flux ratio to the weighted average 
of the K-band flux ratio measurements of all the other epochs 
and fit only the separation, as there is no evidence of statisti- 
cally significant photometric variability over the course of our 
observations. Uncertainties for all parameters are determined 
by fitting each of the individual, stack- averaged images that 
contribute to the final image by the same procedure, and then 
taking the RMS of those values with respect to the average 
value; for the 2002 epoch we also account for the uncertainty 
in the weighted- average flux ratio, which is taken to be the 
RMS of the individual K-band flux ratio measurements at all 
other epochs. 

2.2. Adaptive Optics Data 

TWA 5Aab was also observed using the Keck II 10 m tele- 
scope with the AO system (Wizinowich et al. 2000) and the 
facility near-infrared camera, NIRC 2 (K. Matthews et al., in 
prep), on 2005 February 16 and again on 2005 December 12. 
For these measurements, we used observations of the Galac- 
tic Center to establish that NIRC 2's narrow camera (which 
we used) has a plate scale of 9.961 ± 0.007 mas pixel -1 and 
columns that are at a PA of -0.015 ± 0.134° relative to North 
(J. Lu et al., in prep). In February, five images, each of 0.2 
second exposure time and 30 coadds, were taken through the 
Fell narrow band pass filter (X = 1.65 /im, AA = 0.03 /mi). 
In December, six, six, and three images of 0.181 second ex- 
posure time and 50 coadds were taken through the J (\ = 
1.248 /im, AA = 0.163 /im), H (X = 1.633 /im, AA = 0.296 
/im), and K-prime (X = 2.124 /im, AA = 0.351 /im) band 
pass filters, respectively. These images were processed using 
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Fig. 2.— NIRC 2 AO images of TWA 5A taken on 2005 February 16 
and 2005 December 12. In all images, north is up and east is to the left. 
Component TWA 5 Aa is to the southwest and component TWA 5 Ab is to the 
northeast. 

the same standard data reduction techniques listed above, and 
then shifted and combined to produce a final image. 

Astrometry and flux ratios were then obtained using the 
IDL package StarFinder (Diolaiti et al. 2000). The wide 
brown dwarf companion to the system, TWA 5B, was visi- 
ble in all AO images taken and thus was used as the empirical 
point- spread function required by the StarFinder fitting algo- 
rithm. Both components of the close binary were successfully 
fit by StarFinder in all cases. Errors were calculated by fitting 
the components in all individual images that contributed to the 
combined images and finding the RMS of the values derived 
therein. 

3. RESULTS 

Figure 1 shows the calibrated power spectra and our result- 
ing fits for all five speckle measurements of TWA 5Aab and 
Figure 2 displays our 2005 AO images. We supplement our 
observations with the following measurements that also spa- 
tially resolved TWA 5Aab: the original discovery measure- 
ment by Macintosh et al. (2001), and another taken two days 
later by Brandeker et al. (2003). Table 1 lists all separation, 
position angle, and flux ratio measurements that are used in 
this study. 

Over the six years that the components of TWA 5 A have 
been spatially resolved, the binary has undergone a full or- 
bit (see Figure 3), allowing an accurate estimate of its orbital 
parameters. We calculate an orbital solution for TWA 5 A us- 
ing the Thiele-Innes method (e.g., Hilditch 2001), minimizing 
the x 2 between the model and the measurements, which are 
converted from angular separation and position angle to right 
ascension and declination. Our model incorporates the fol- 
lowing 7 standard orbital elements: P (period), A (semi-major 
axis), e (eccentricity), i (inclination), T (time of periapse pas- 
sage), Q (longitude of the ascending node), and uo (argument 
of pericenter). With nine two-dimensional astrometric mea- 
surements, there are eleven degrees of freedom in our fit. The 
best fitting model produced a x 2 of 8.91 with 11 degrees of 
freedom. Furthermore, of the nine data points used in the fit, 
seven are within la and two are within 2a of the model, sug- 
gesting that our fit is good. The la uncertainties in the model 
parameters are estimated by changing the values of x 2 by one 
(Bevington & Robinson 1992). Table 2 lists the best-fit orbital 
parameters and their uncertainties and Figure 3 shows our so- 
lution for the projected orbit of TWA 5 A. This astrometric so- 
lution yields a mass 0.71 ± (0.14 ± 0.19) M©, where the two 
sources of uncertainty stem from our orbital solution and the 
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TABLE 2 
Orbital Parameters of TWA 5 A 



P (years) 5.94 ± 0.09 

A (") 0"066 ± 0.005 

T (years) 2004.34 ± 0.09 
e 0.78 ± 0.05 

i (degrees) 97.4 ±1.1 

Q (degrees) 37.4 ± 0.9 

uj (degrees) 255 ± 3 



Note. — We note that the value of Q is actually subject to a 180 degree ambiguity 
without three dimensional velocity information. 
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Fig. 3.— Orbital solution for TWA 5Aa and TWA 5Ab. Filled circles 
represent the data taken in this study, while the open circles represent data 
taken from the literature (both in 2000). The star at [0,0] marks the position 
of TWA 5Aa; the lines from the data points to the ellipse indicate where the 
fit believes the point should lie on the orbit. The dashed line through the 
center represents the line of nodes and the 'X' marks the location of closest 
approach. The parameters for this orbit are given in Table 2. 

from the 9% uncertainty in the kinematic distance estimate to 
TWA 5A (Mamajek 2005). 

4. DISCUSSION AND CONCLUSIONS 

With a well-determined astrometric orbit for TWA 5A, it 
is possible to begin to compare dynamical estimates of mass 
with those inferred from theoretical PMS tracks. By compar- 
ing the quantity M/ (D/44 pc) 3 , we preserve the high precision 
of the astrometric data in the analysis. Here we investigate 
the models by Baraffe et al. 1998 (a = 1.0), D'Antona & 
Mazzitelli 1997, Palla & Stahler 1999, Siess et al. 2000, and 
Swensonetal. 1994. 

Model estimates of mass and age require both effective tem- 
perature and bolometric luminosity as inputs. Effective tem- 
peratures are estimated from the unresolved spectral type for 
TWA 5 A and component flux ratios. The photometric analy- 
sis performed here and elsewhere (see Table 1) shows that the 
two components of TWA 5 A have nearly equal brightness out 
to 1.6 fim. We therefore assume that the components have the 
same spectral type and assign it to be Ml. 5 ± 0.5, the spec- 
tral type found from a spatially unresolved spectrum (Webb 
et al. 1999). This is consistent with the J-H colors for each 
component, which are calculated by combining the 2MASS 
unresolved magnitudes with the flux ratios measured here (J- 



H[M1.5] = 0.67; Leggett et al. 1992). We estimate the tem- 
perature of these components using a conversion of spectral 
type to effective temperature given in Luhman et al. (2003). 
These temperatures are intermediate between dwarf and giant 
stars, thus making them appropriate for PMS stars like TWA 
5A. Bolometric luminosities are estimated using our H band 
magnitudes and the corresponding H band bolometric correc- 
tions for PMS stars (Luhman 2005, private communication), 
along with the kinematic distance for TWA 5A of 44 ± 4 pc 
(Mamajek 2005). These input values are given in Table 3. 

The uncertainties in M/ (D/44 pc) 3 and age from each 
model are estimated through Monte Carlo simulations. Our 
Monte Carlo simulation was run with 10 5 points sampled 
from random, independent gaussian distributions of H-band 
flux densities and flux density ratios, temperatures for each 
component, and distance. This allowed us to correctly ac- 
count for the correlation between input values. In particular, 
the bolometric correction used for determining the luminosity 
stems from the effective temperature. Additionally, the lu- 
minosities of the primary and the secondary are correlated, as 
they are calculated from the same parameters, namely the total 
flux and the flux ratios. Each of of these values are converted 
into sets of temperatures and bolometric luminosities, which 
are is fed into the models to produce estimates of masses and 
ages for each component. Using the distance assumed for 
each run, we convert the total mass of both components to 
M/ (D/44 pc) 3 . We emphasize that though the predictions of 
the models are clearly not independent of the distance uncer- 
tainty, by using this quantity as opposed to the total mass, 
we can leave the distance uncertainty out of our dynamical 
measurement. Thus, we then generate contours that correctly 
represent the la uncertainties on the masses and ages derived 
from the model tracks. 

Figure 4 shows the dynamical mass normalized by the 
distance-cubed and with it the estimated age of the TW Hy- 
drae association (10 ± 2 Myr, based on comparison with the (3 
Pictoris moving group, Zuckerman & Song 2004), compared 
with the model predictions as determined by the Monte Carlo 
simulations. While the input values are generally treated as 
independent variables, the Monte Carlo simulations demon- 
strate that it is important to consider the existing correlations 
between input values, as the resulting uncertainties are oth- 
erwise significantly underestimated (by up to a factor of five 
depending upon the model). With this proper treatment, the 
dynamical mass and age estimates are within 2a of all five 
tracks. 

While it is reassuring that current models agree with these 
measurements, we are unable to distinguish between the 
model tracks. Future improvements in the precision of both 
the distance and the temperature would allow for a more il- 
luminating comparison between the theoretical models. The 
temperature uncertainties on the components of TWA 5A 
could be substantially reduced with spatially resolved spectral 
types, which should drive the uncertainty in the temperature 
down to ± one quarter of a spectral subclass. Uncertainties 
in the distance could also be substantially reduced. Typical 
uncertainties in new parallax measurements are on the order 
of a few milliarcseconds, with improvements on these results 
promised in the near future (Vrba et al. 2004). A separate 
method to determine the distance is the addition of spatially 
resolved radial velocities, which would allow an independent 
distance estimate from the orbital parameters. In principle, a 
factor of two improvement in the distance uncertainty could 
be achieved. These radial velocity measurements would also 
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TABLE 3 
TWA 5 System Photometry 



Component 


Sp. Ty. 


m K 


m H 


my 


J-H 


H-K 


Log[Temp] (K) 


Log[Lum] (L ) 


TWA 5Aa 
TWA 5Ab 


M1.5 ±0.5 
[Ml. 5 ±0.5] 


7.39 ± 0.04 
7.62 ± 0.08 


7.69 ± 0.04 
7.79 ± 0.05 


8.40 ± 0.07 
8.45 ±0.15 


0.71 ± 0.08 
0.66 ±0.16 


0.30 ± 0.06 
0.17 ±0.09 


3.56 ± 0.01 
3.56 ± 0.01 


-0.80 ± 0.08 
-0.84 ± 0.08 



Note. — The apparent K, H, and J band magnitudes for TWA5Aa and Ab are computed using 2MASS measurements of the combined magnitudes and the flux ratio measurements 
given in Table 1 . For the H and J band, we use the most recent AO flux ratio measurements taken with NIRC 2 in 2006 December. For the K band, we use the average flux ratio 
from all speckle measurements (excluding the unresolved measurement in 2002). The spectral type for TWA 5Ab is assumed to be roughly the same as TWA 5Aa, given that the two 
components have nearly equal brightness. Temperatures are estimated from the scalings given in Luhpan et al.^20T^3, while^ol^me^jic luminosities gtg^al|gl^ed usiiy^the kij^matic 



distance and H 



band bolometric corrections (Luhman 2005, private communication) . 
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AD'Antona &c Mazzitelli 1997 
T Palla & Stahler 1999 
^Siess et al. 2000 
■ Swenson et al. 1994 



ever, no mid-infrared excess has been detected in unresolved 
measurements. Given the tightness of the TWA 5 A binary, 
any disk material is likely to be localized to a very radially 
thin reservoir of material near the dust sublimation radius. In 
either case, disk or low mass companion, the cause of the in- 
frared excess should not have a large impact on the track com- 
parison, as the comparison was performed at H-band, where 
the luminosities do not appear to be significantly effected by 
this excess. 

In summary, our solution to the orbit of this system and the 
subsequent determination of its mass shows that these young, 
nearby associations of stars are excellent laboratories for the 
study of low-mass star formation. There are likely other sys- 
tems in TW Hydrae with similar close companions that will 
yield more mass estimates within a short time, much like we 
have seen here. Thus monitoring of these systems will greatly 
aid in the constraint of PMS mass tracks in the near future. 



Fig. 4. — Mass and age for TWA5Aa+b calculated from the dynamical so- 
lution and predicted from each of the five theoretical models considered. The 
age shown here is the age predicted for the primary, but the age predictions 
for both components are consistent with each other. The masses are divided 
by the distance cubed. Though the predictions of the models do include the 
distance uncertainty, by using the quantity M/ (D/44 pc) 3 , we insure that only 
the model contours are affected by distance uncertainty, avoiding correlation 
between the models and our dynamical measurement. 

improve the total mass estimate, as would additional astro- 
metric measurements, and would also eventually allow for 
the determination of individual masses of each of the com- 
ponents. However a potential complication is the possible ex- 
istence of an additional component in the system (Torres et al. 
2003). Another possible method of determining the individ- 
ual masses of the components of TWA 5A is to use absolute 
astrometry with respect to TWA 5B, as was done recently for 
T Tau S (Duchene et al. 2006). If we just take the improve- 
ments expected from future parallax measurements, spatially 
resolved spectra, and additional astrometry, we expect to be 
able to distinguish between the tracks at the ~ 3 a level. 

Future spatially resolved spectral types would also address 
the marginal inconsistency between the J-H and H-K colors 
for the primary. Currently, we have assumed that this arises 
from a small K-band excess (2a). Torres et al. (2003) suggest 
that an additional component may be present in the system. If 
this component were particularly low mass, it could give rise 
to the apparent K band excess of the primary. Alternatively, 
the infrared excess could arise from circumstellar material. 
Mohanty et al. (2003) report the detection of strong Ha emis- 
sion in an unresolved, high resolution spectrum of TWA 5A, 
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APPENDIX 
NIRC PLATE SCALE AND ORIENTATION 

NIRC's pixel scale and orientation are calibrated relative to NIRC 2 using observations of the Galactic center (see Table 5). 
As discussed in section §2.2, NIRC 2's absolute offset with respect to north has been accurately calibrated, so that by calibrating 
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TABLE 4 

Absolute NIRC Position Angle Offsets and Uncertainties 



Epoch(s) 


No. of 


Absolute 




Measurements 


PA (deg) 


1998 April - 1998 August 


4 


-0.40 ±0.135 


1998 October - 2002 July 


11 


-0.884 ±0.143 


2003 April - 2003 September 


3 


0.761 ±0.135 


2004 April - 2005 May 


4 


-0.728 ±0.196 



NIRC with respect to NIRC 2, we can use a simple coordinate transform to get its absolute orientation as well. While NIRC's 
pixel scale has been stable, known engineering adjustments have introducted slight rotations in the camera over time. Four out of 
five of TWA5A speckle measurements were taken at a time with NIRC's orientation is well-characterized by the Galactic center 
experiment. The 2003 measurement, however, was taken during a period of multiple engineering adjustments and no Galactic 
center data. We therefore bound its orientation by the measurements taken of the Galactic Center just before and after it, which 
leads to an absolute offset of 0.032 ± 0.719°. For all of our observations, we use a constant plate scale of 20.45 ± 0.03 mas/pixel. 
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